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Abstract 

In this review paper we consider the polynuclear growth (PNG) 
model in one spatial dimension and its relation to random matrix 
ensembles. For curved and flat growth the scaling functions of the 
surface fluctuations coincide with limit distribution functions coming 
from certain Gaussian ensembles of random matrices. This connection 
can be explained via point processes associated to the PNG model and 
the random matrices ensemble by an extension to the multilayer PNG 
and multi-matrix models, respectively. We also explain other mod- 
els which are equivalent to the PNG model: directed polymers, the 
longest increasing subsequence problem, Young tableaux, a directed 
percolation model, kink-antikink gas, and Hammersley process. 

1 Introduction 

In this paper we consider a stochastic growth model, the polynuclear growth 
(PNG) model, on a one- dimensional substrate. The relevant features of the 
dynamics of this model stochastic local growth rule with a smoothing 

mechanism [54^ . The latter prevents the formation of large spikes. As a con- 
sequence the surface, on a macroscopic scale, follows a deterministic growth 
rule: it has a limit shape. Nevertheless, on a mesoscopic scale the surface 
is still rough. This roughness is the observable we are mainly interested in. 
The PNG model belongs to the KPZ universality class. The KPZ model 
was introduced by Kardar, Parisi, and Zhang in a seminal paper [3T] where 
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they described random surface growth by a non-linear stochastic differential 
equation. For details on the universality we refer to Prahofer's thesis |55j . 
Chapters 2 and 3, see also pHj- In one dimension one can determine the dy- 
namical exponent, z = 3/2, by a scaling argument or renormalization group 
methods, see for example the books [SH EH]. As a consequence, for large 
growth time t, the fluctuations of the surface height scale as t 1 / 3 and the 
spatial correlation length is of order t 2 ^ 3 . These exponents should hold for 
all growth models in the KPZ universality class. It is commonly expected 
that not only exponents but also scaling functions and limiting distributions 
are universal. To study these more detailed informations on KPZ growth, 
one looks for simplified but still solvable models in the KPZ class. Hence the 
study of the PNG model. 

A main step forward in understanding KPZ growth and suggesting a 
deep connection to random matrices was achieved by Johansson jHEj- He 
considered a discrete growth model on a one-dimensional substrate, which 
can be interpreted, among others, as a kind of first-passage site percolation 
model JHl- For specific initial conditions, he obtained the surprising link 
between the shape fluctuations of the percolated region and the GUE Tracy- 
Widom distribution, F 2 , which was first introduced in the random matrix 
context [OH]- Shortly before that, the same distribution appeared in a work 
by Baik, Deift, and Johansson 6 on the problem of the longest increasing 
subsequence in a random permutation. On the other hand there was the 
PNG model, a model considered by physicists since the 60's. It was known 
that it belongs to the KPZ class. Then Prahofer and Spohn jHE] noticed the 
mapping between (a version of) the PNG model and the increasing subse- 
quences problem. Hence the connection between KPZ and random matrices. 

Here we mainly focus on two geometries of the PNG model: curved growth 
and flat growth. The first generates a droplet-like profile, hence called PNG 
droplet. The macroscopic profile of the second is flat, thus flat PNG. For these 
two geometries more refined information is available. In the limit of large 
growth time, properly rescaled, the height fluctuations of the PNG droplet are 
described by the F 2 distribution [HHl IH] - F2 is the limiting distribution of the 
properly rescaled largest eigenvalue of the Gaussian Unitary Ensemble (GUE) 
of random matrices. Moreover, the limiting process describing the surface 
height of the PNG droplet has been identified as the Airy process |SZ|, which 
also arises in the multi-matrix extension of GUE. It describes the evolution 
of the largest eigenvalue of Dyson's Brownian motion. 

For the second geometry considered, the flat PNG, the picture is not 
yet complete. It is known [SHI E] that the scaling function of the height 
fluctuation is the GOE Tracy- Widom distribution, F\ [ZD]- This distribution 
arises in the Gaussian Orthogonal Ensemble (GOE) of random matrices. F\ 
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is the limiting distribution of the largest eigenvalue of GOE. It is conjectured 
that it should correspond to the evolution of the largest eigenvalue of Dyson's 
Brownian motion (for orthogonal matrices). A result going in the direction 
of this conjecture is obtained by Ferrari (2D]- The problem remains open 
for the GOE multi-matrix model. Very recently, Sasamoto [BH! obtained a 
process in the context of the totally asymmetric exclusion process, which, by 
universality, should describe also the surface height of the flat PNG. 

The connection between the PNG model and the Gaussian ensembles of 
random matrices can be understood via point processes. Although random 
matrices are not directly related to the PNG model, it turns out that both 
can be described by point processes with the same mathematical structure. 
Under an appropriate scaling one obtains the same limit point processes when 
the growth time, resp. matrix dimension, tends to infinity. 

The paper is organized as follows. In Section El we introduce the PNG 
model and report known results. In Section El further equivalent models are 
described. In Section0]we explain the extension to the PNG multilayer model 
and the related point processes. The link between line ensembles and (real- 
valued) Young tableaux is also discussed. In Section El we introduce Gaussian 
ensembles of random matrices, the point processes of their eigenvalues, and 
the extension to multi- matrix models. Section |H1 is devoted to the discussion 
of the connection between random matrices and the PNG model as well as 
Young tableaux. 

2 The polynuclear growth model 

Let us describe the polynuclear growth (PNG) model in continuous time 
in 1 + 1 dimension. It is a growth model on a one-dimensional substrate. 
The surface at time t is described by an integer-valued height function x i— > 
h(x,t) G It, see Figure [TJ To be precise, at the discontinuity points of h, the 
height function has upper limits, i.e., {x G H\h(x, t) > k} is a closed set for 
all fceZ. For fixed time t G R, consider the height profile x i— > h(x,t). The 
height h, as x increases, has jumps of height one at the discontinuity points, 
called up-step if h increases and down-step if h decreases. Finally, if at time 
t there is a spike, i.e., a pair of up- and down-steps at the same position x, 
then we call it a nucleation event at (x, t). 

The PNG dynamics has a deterministic and a stochastic part: 

(a) Deterministic part: the up-steps move to the left with unit speed and 
the down-steps to the right with unit speed. When a pair of up- and 
down-steps collide, they disappear. 
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Figure 1: The PNG dynamics. The up- (down-)steps move to the left (right) 
with unit speed. The big arrow represents a nucleation. 



(b) Stochastic part: the nucleation events form a locally finite point process 
in space-time (usually a Poisson process). Once a pair of up-down steps 
is created, it immediately follows the deterministic dynamics. 

The stochastic part of the dynamics produces the roughness of the surface. 
This is counterbalanced by the smoothing due to the deterministic part. 

Typically one considers flat initial conditions, h(x, 0) = for all x G 
R, with the nucleation events given by a Poisson process (not necessarily 
with uniform intensity). By varying the space-time intensity g(x,t) of the 
Poisson process, different geometries can be obtained. Below we consider 
two cases of particular interest. For the PNG droplet the nucleations occur 
with constant intensity in the region spreading with unit speed from the 
origin in both directions. For the flat PNG the nucleations have constant 
intensity everywhere, thus the surface is statistically translation-invariant. A 
visualization of these two geometries is provided at ^H] as a Java applet 

At this point it is convenient to introduce some vocabulary taken over 
from special relativity. The steps move with speed of light ±c, c = 1. So 
their trajectories in space-time have slope ±1, called light-like. The forward 
light cone of a point (x, t) is the set of points {(a/, t') \ \x — x'\ < t! — t}, and 
the backward light cone of (x,t) is {(x',t')\ \x — x'\ < t — t'}. With light 
cone we denote the union of the forward and the backward cone. A path is 
called time-like if it is included in the light cone of each of its points. A path 
is called space-like if the light cone attached to each of its points does not 
contain any other point of the path. 
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2.1 The PNG droplet 

The PNG droplet is obtained from a flat initial height profile, h(x, 0) = 
for all x G R, if the density of Poisson points is constant (here we choose 
g — 2) in the forward light cone of the origin and zero outside, i.e., for 
(x, t) e R x R+ 





if \x\ 


{o 


if x 



<?(*.*) H n ;*L (2.1) 



For large growth time t the typical shape of the PNG droplet is a half circle, 
see Figure |2l More precisely, it converges to the limit shape given by 

lim t _1 /i(rt,t) = 2Vl - r 2 , re [-1,1], (2.2) 

in probability. 

To see the roughness of the surface one has to look at a mesoscopic scale 
around the macroscopic shape 2ty/l — (x/t) 2 , x G [—£,£]. A first natural 
question concerns the scale of fluctuations and their limit behavior. The 
result is that the vertical fluctuations live on a t 1 ^ 3 scale with limiting distri- 
bution function 

lim P(h(0,t) < 2t + t 1/3 s) = F 2 (s), (2.3) 

t— >oo 

where F 2 is the GUE Tracy- Widom distribution function The conver- 
gence is in distribution as well as for all finite moments. ()2.3j) was obtained 
in [HE] by mapping the PNG droplet to the Poissonized version of the prob- 
lem of longest increasing subsequences 16]. Similarly if one looks away from 
x = one has, for any fixed r G (—1,1), 



lim P(h(rt,t) < 2tVT^ + t 1/3 (l -r 2 ) 1/6 s) = F 2 {s). (2.4) 
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The second interesting question concerns the spatial height correlations. 
The correlation length scales as t 2//3 for large growth time t. Therefore one 
defines the rescaled surface height as 



£ -> h? SC (0 = t- 1/3 (M^ 2/3 ,i) - 2t V / l-^- 2/3 )- (2-5) 
In jHZj it is shown that, in the sense of finite dimensional distributions, 

lim hr c (£)=A(0, (2.6) 



t- 



where A is the Airy process. The definition and properties of the Airy process 
are given in Section 14.31 This process arises also in the multi-matrix model 
for GUE random matrices as we will explain in Section 1531 

More recently Borodin and Olshanski showed pH] that the Airy process 
describes the space-time correlations not only for the space-time cut with 
constant t, but also along any space-like (and light-like) path in the droplet 
geometry. They prove convergence of finite-dimensional distributions in the 
language of Young diagrams. The link with Young diagrams is explained in 
Section f4.21 For each point (u, v) G R+, they consider the random Young 
diagram Y(u,v) obtained by the RSK correspondence. Then for each space- 
like path in R+ a Markov chain is constructed, which describes the evolution 
of the Young diagram Y. The case u+v = t is the one of the PNG droplet [5T] . 
The case uv = constant corresponds to the terrace-ledge-kink (TLK) model 
which can be regarded as model for the facet boundary of a crystal in thermal 
equilibrium [231 121]- For time-like paths no result is known. The major 
difficulty seems to lie in the lack of a Markov property. 



2.2 Flat PNG 

A flat initial condition, h(x, 0) = for x G R, and constant density of 
Poisson points in R x R + (as before we choose g = 2) generates the flat PNG 
geometry, see Figure El Since no other constraint is imposed the surface 
height h(x,t) is statistically translation- invariant, thus we focus on x = 0. 

The fluctuations live on the same scale as for the PNG droplet, namely 
t 1 / 3 . However, the limiting distribution is different as shown in [HE], namely 

lim P(/i(0,t) < 2t + i 1/3 2~ 2/3 s) = iq(s), (2.7) 

t— >oo 

with F\ the GOE Tracy- Widom distribution function [SH- This result fol- 
lows from a related problem on longest increasing subsequences [H]. The 
convergence is in distribution and for all moments as in ()2.3|) . 
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Figure 3: A sample of the flat PNG. 



In [20] a first step towards understanding the spatial correlation is made, 
as we will explain in Sections 14. 31 and IHTT1 In a recent preprint, Sasamoto 
obtained the spatial correlation for a closely related model. By universality 
this new process should be the analogue of the Airy process for flat PNG. 

2.3 A geometrical point of view for the PNG model 

Let us illustrate the space-time picture of the PNG model, which is the 
starting point for the link to other equivalent models described in the next 
section. 

First we consider the case of flat initial conditions, h(x, 0) = for all 
x G R. For a given realization of Poisson points, we construct h(x, t) for 
t G [0, T], T > fixed, as follows. One starts plotting the nucleation events 
in space-time. Then, increasingly in time, one draws the trajectories of the 
up- and down-steps in space-time. These are light-like paths. When two of 
these paths meet, as t increases, they stop. This reflects the disappearing 
of the corresponding up- and down-step. In this way one divides space-time 
into regions bounded by piecewise straight lines with slopes ±1, see Figure EJ 
The height h(x, t) is constant in each region and is given by the number of 
lines crossed by any time-light path from (x, 0) to (x, t). For any given t, 
since the Poisson process is a.s. locally finite, so is the number of steps of 
x I— > h(x,t), thus h is a.s. locally bounded. 

In the case that the initial surface profile is not flat, the surface height 
at some later time t is obtained in a similar way. The only difference is 
the following. To the lines generated by the Poisson points we need to add 
additional lines starting from the t — axis with slope —1, resp. +1, if 
initially at x there is an up-step, resp. a down-step. Then, the number 
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of lines crossed along any time-like paths from (x, 0) to (x, t) is the height 
difference h(x,t) — h(x,0). 

2.4 Other geometries studied for the PNG model 

As explained above, the one-point distribution function of the surface height 
for both geometries, the PNG droplet and flat PNG, was obtained by iden- 
tifying the surface height with the length of a longest directed polymer. The 
directed polymer can be mapped to a longest increasing subsequence in a 
random permutation (without /with involution), for which Baik and Rains 
have obtained the asymptotics |H1 Ej by using Riemann-Hilbert techniques. 
For the PNG droplet, the joint-distribution of the height profile is obtained 
in [2Zj using a multilayer generalization of the PNG as will be explained 
in Section 0] The multilayer technique has been motivated by the work of 
Johansson on the Aztec diamond I4()j. 

Sasamoto and Imamura studied the (discrete) half-droplet PNG geome- 
try, which consists in allowing nucleations only in the region x G [0,t] |33j . 
They prove that the rescaled height is GUE distributed away from x = and 
that there is a transition to GSE at x — 0. If extra nucleations are added at 
the origin with intensity 7 > 0, the distribution above x = has a transition 
at 7 = 1. For 7 < 1 it is still GSE, for 7 = 1 it is GOE distributed, and 
for 7 > 1 the fluctuations become Gaussian, because the contribution of the 
nucleations at the origin dominates. The asymptotic one-point distributions 
at the origin follow from jUlE], too. 

A modification of the PNG droplet consists in adding sources at both 
extremities of the droplet. Extra nucleations with fixed line density a + and 
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a_ are independently added along the boundary of the forward light cone of 
the origin, i.e., in (x,t) such that \x\ = t. This model was used in jSH] to 
describe stationary PNG growth. Baik and Rains [Z| obtained detailed results 
for the asymptotic distributions which we describe briefly. For a± small, the 
effects coming from the edges are small and the fluctuations are still GUE 
distributed. On the other hand, if a + > 1 or a_ > 1, the boundary effects are 
dominant and the fluctuations become Gaussian. The cases where a + — 1 
and/or a_ = 1 are also studied and other statistics arise. Of particular 
interest is when = 1 for 1 — a± — £>(r 1/3 ), in which case the PNG 

growth is stationary and has a flat limit shape. Stationary PNG together 
with flat PNG and PNG droplet are the three most interesting situations in 
the context of surface growth. 

Johansson in j^H] pointed out the connection between random matrices 
and the shape fluctuations in a discrete model of directed polymers, equiva- 
lent to the discrete PNG model. A class of corner growth models in discrete 
space-time (similar to the PNG droplet) is analyzed in [301133]. The discrete 
versions of the different geometries for the PNG model discussed above are 
studied in a series of papers jH3 HUH EH EH] • 

3 Equivalent models 

In this section we discuss the mapping between models which are equivalent 
to the PNG model. We start explaining the directed polymers on Poisson 
points. The link to the PNG was used in jHEj to obtain the first results on the 
height fluctuations of the PNG model. Then we continue with the longest in- 
creasing subsequence problem, Young tableaux, a directed percolation model, 
the kink-antikink gas, and the Hammersley process. 

3.1 Directed polymers on Poisson points 

Let us define a partial ordering in R 2 , -<, as follows. For y, z G R 2 we say 
that y -< z if both coordinates of y are less or equal than those of z. Consider 
a Poisson process with intensity 1 in R 2 . For a given realization of Poisson 
points, a directed polymer on Poisson points starting at (0, 0) and ending at 
(t, t) is a piecewise linear path 7 connecting (0, 0) -< q% -< . . . -< g^ 7 ) -< (t, t), 
where qi are Poisson points. The length ^(7) of the directed polymer 7 is the 
number of Poisson points visited by 7. The basic observable of interest is the 
maximal length of the directed polymers, 

L(t) = maxZ(7). (3.1) 

7 
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Figure 5: Height and directed polymers for the droplet geometry 



This is called point-to-point setting because both initial and final points are 
fixed. 

A modification of the problem consists in considering the set of directed 
polymers starting from (0,0) and ending in the segment XJ t = {(y, z) G 
IR+ly + z = 2t}. This is called the point-to-line problem and the maximal 
length is denoted by L^it). 

The link with the PNG model is apparent once we use the graphical point 
of view explained in Section T2.3I In fact, for the PNG model, h(0,t) equals 
the number of lines (up- and down-steps trajectories) crossed by any light- 
like paths from (0,0) to (0, t). In particular, one considers the paths which 
cross them at the nucleation points and consist in straight segments between 
these points. These are the directed polymers introduced above, up to a 7r/4 
rotation, see Figure El Because of the 7r/4 rotation, if the density of Poisson 
points in the PNG model and the directed polymer is the same, then h(0,t) 
equals L(t /y/2). To have a nicer formula, we set the density of Poisson points 
to 2 for the PNG model and to 1 for the directed polymers. This implies 
that L(t) = h(0,t) of the PNG droplet and L e (t) = h(0,t) of the flat PNG, 
both in distribution. Thus 



lim P(L(t) <2t + st 1/3 ) = F 2 (s), lim P(L t (t) < 2t + s2~ 2/3 t 1/3 ) = F x (s). 

t^oo t^oo 

(3.2) 



3.2 Longest increasing subsequences 

Let Sn denote the permutation group of the set {1, . . . , iV}. For each per- 
mutation a G Sn, the sequence (cx(l), . . . , cr(N)) has an increasing subse- 
quence of length k, (a(ni), . . . ,a(rik)), if <j(m) < o-(n 2 ) < . . . < cr(nfc) and 
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rii < n 2 < . . . < n k . Denote by L N (a) the length of the longest increasing 
subsequences for the permutation a. The problem of finding the asymptotic 
law of Ln for a uniform distribution on Sn is also called Ulam's problem 
(1961) [25] • For a review around this problem, see j3j. Baik, Deift and Jo- 
hansson in the seminal paper |Hj determined the fluctuation law of L^. They 
prove 

lim P(L N < 2y/N + sN 1/6 ) = F 2 (s), (3.3) 

N^oo 

where F 2 is the GUE Tracy- Widom distribution function. Compare ()3.3)1 
with (|3.2jl . the role of t is taken over by \/N. 

()3.3|) is obtained using the Poissonized version of the problem, which is 
the approach of the problem used by Hammersley [32] ■ Instead of fixing 
the length of the permutations to N, one considers the set of permutations 
S = U n >o<S n and assigns the probability e~ N k N jk\ that a permutation is 
in Sk- Baik et al. first prove that ()3.3j) holds for this problem, and then 
obtain the result via a de-Poissonization method, consisting in bounding 
from above and below the distribution of in terms of the Poissonized one. 
In a statistical physics language, the problem with fixed iV corresponds to 
the canonical ensemble, the one with Poisson distributed length to the grand 
canonical ensemble, and the N — > oo limit to the thermodynamical limit. It 
is not surprising that the equivalence of ensembles holds for the observable 
Ljy. In fact, in the grand canonical ensemble the typical value of the number 
of Poisson points is 0{\/N) apart from N. Thus the correction to is of 
order 1, which is vanishing small compared to iV 1 / 6 as iV — > oo. 

The problem of the longest increasing subsequences in Sn is equivalent 
to the problem of finding the longest directed polymer from (0, 0) to (t, t) 
when N points are distributed uniformly in the square [0,t] 2 . The directed 
polymers on Poisson points is the Poissonized version. In fact, consider a 
configuration of iV points in the square [0,t] 2 . The length of the longest 
directed polymer depends only on the order of their projections along both 
axis. Without changing this order we can put them on {1, ... , iV} 2 . Thus to 
each directed polymer there corresponds an increasing subsequence with the 
same length. Hence the length of the longest directed polymer equals the 
one of the longest increasing subsequence. See Figure EH for an example. 

3.3 Young tableaux 

Let A = (Ai, A 2 , • • • , Afc) be a partition of the integer N, i.e., satisfying 
Ai > A 2 > . . . Afc > 1 and Yli=i Ai = A^. A is represented by a diagram with 
k rows and with Aj cells for row i, called Young diagram. A (standard) 
Young tableau of shape A is a Young diagram, where the cells are filled 
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Figure 6: Longest increasing subsequences and directed polymers. (2, 3, 5) 
and (2,3,4) are the two longest increasing subsequences. The directed poly- 
mers of maximal length passes through the points labelled by (2, 3, 5) and 
(2,3,4). 



by the numbers 1, 2, . . . , N, increasingly in each row and column. The 
Robinson-Schensted correspondence is a bijection between permutations 
cr e Sn and pairs of Young tableaux (V(a), Q(cr)) with N cells and the 
same shape. The algorithm leading to (V(cr), Q(<j)) is the following 

One starts with a pair of empty tableaux V and Q and fills the 
cells as follows: 

P-tableau: for i from 1 to N, the number a(i) is always placed 
via row-bumping in the first row, i.e., 

a) if a(i) is larger than all numbers in the first row of the 
P-tableau, then append to the right of them, 

b) otherwise put it at the place of the smallest entry in the first 
row of V, which is larger than u{i). 

In case b), the entry which was replaced is now placed via 
row-bumping in the second row. If an entry is replaced in the 
second row, then it is placed via row-bumping in the third row, 
and so on. 

Q-tableau: At each step in the generation of the P-tableau, its 
shape is enlarged by one cell. For each i from 1 to N, put the 
number i at the position where a new cell appeared at step i in 
the P-tableau. 
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5 


a(i) 


2 


3 


1 


5 


4 


V 


2 


2 3 


1 3 
2 


1 3 5 
2 


1 3 4 

2 5 


Q 


1 


1 2 


1 2 
3 


1 2 4 
3 


1 2 4 

3 5 



Table 1: Construction of the Young tableaux V(a) and Q(cr) for the permu- 
tation a = (2,3,1,5,4). 



As an illustration in Tabled we show the construction of the Young tableaux 
for the permutation a = (2,3,1,5,4) of Figure E3 By construction, the 
Young tableaux V(a) and Q(er) have the same shape. Given a permutation 
o G <S>tv, one places the points on N 2 with coordinates (i, cr(i)) and draws the 
lines as in Figure H3 Notice that the first row of V(a) contains precisely the 
positions of the horizontal lines at abscissa N + 1/2, likewise the first row of 
Q(er) contains the positions of the vertical lines at ordinate N + 1/2. This 
property holds for all permutations which means that [I] 

L N (a) = Ai(<r). (3.4) 

Consequently, a way to determine the asymptotic behavior of L N is by 
analyzing the length of the first row of Young tableaux [T^J El] . The measure 
on the set of partitions of {1, ... , N}, Y N , induced by the uniform measure on 
Sn via the RS correspondence is the Plancherel measure Pbv; let d\ denote 
the number of Young tableaux of shape A, then 

PLv(A) = dx A G Y N . (3.5) 

The lengths A2(er), A3 (a), . . . also have an interpretation in terms of di- 
rected polymers. This is discussed in Section l4~2l 

3.4 A directed percolation model 

Motivated by the discrete anisotropic directed polymers model of Rajesh 
and Dhar jnHl, one can reformulate the PNG model starting from a directed 
polymer picture as follows. We take a stack of horizontal planes R 2 xN C R 3 . 
Neighboring planes are connected by vertical bonds placed randomly. One 
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introduces the percolation cluster P C R 2 x N as follows. In each plane there 
is perfect directed percolation, i.e., if (x, y,k) G P then {(x 1 , y', k)\x' > x, y' > 
y} C P. Percolation between adjacent planes occurs through the bonds, i.e., 
if there is a bond from (x, y, k) to (x, y,k + l), then (x, y, k + 1) G P whenever 
(x,y,k) G P. 

If the position of the bonds between planes are given by independent 
Poisson processes of intensity 1, then the cluster spreading from the origin, 
(0,0,0), corresponds to the point-to-point directed polymer. The height of 
the cluster at (x, y) is the largest k such that (x, y, k) G P. Thus the height 
at (x, y) equals, in distribution, the length of the longest directed polymer 
from (0, 0) to (x, y). Similarly, the point-to-line directed polymer corresponds 
to the cluster spreading from the line {(x, y, k)\x + y = 0, k = 0}. Define 
the time axis as {(x, x)\x G Then the PNG dynamics is recovered by 

slicing the percolating cluster perpendicularly to the time axis. 

The percolation cluster arises as the continuum limit of the discrete model 
of Rajesh and Dhar jJJJjj. It appears also in the limit of large alphabets for a 
model of sequence aligning 

3.5 Interacting particle models 
The kink-antikink gas 

Bennet et al. JT] studied an idealized model for solitons in the sine-Gordon 
model. There are two types of solitons, called kinks and antikinks. Kinks 
move to the left with velocity —1, antikinks to the right with velocity 1. 
The solitons are point-like and do not interact, with the exception that, if a 
kink and an antikink collide, they annihilate each other. On the other hand, 
there is a constant uniform rate of production of kink-antikink pairs, which 
immediately move apart with unit speed. In pi] the stationary distribution 
of the kink-antikink gas is studied on a finite ring and in the thermodynamic 
limit. From the description of the PNG model it is clear that the kinks 
and antikinks can be identified with the up- and down-steps of the PNG 
height profile. Thus any result for the PNG model has a translation to 
statements for the kink-antikink gas and vice versa, provided initial and 
boundary conditions match. For example, the height above the origin for 
flat PNG equals the number of kinks and antikinks that passed through the 
origin (with vacuum as initial condition). 
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Hammersley process 

In order to tackle the longest increasing subsequence problem, Hammers- 
ley used a particle system to prove that limjv^oo ~Ei(Ln) /\/N exists, see 
also [31 EE]. There is only one kind of particles on [0,T], which sometimes 
jump to the left a certain distance, but otherwise are at rest. The position 
x = T serves as particle reservoir and initially no particles are in [0, T). If a 
particle jumps, it jumps to the left to a position uniformly chosen in the in- 
terval between the particle and its left neighbor (or, for the left-most particle, 
the origin). The jump rate for each particle is proportional to the distance to 
its left neighbor. Particles do not cross under the dynamics. Moreover, the 
space-time points where the particles land form a Poisson point process in 
[0, T) x R + with uniform density. It is easy to see that for a given realization 
of Poisson points, the space-time trajectories of the particles correspond, up 
to a 7r/4 rotation, to the light-like lines generated by the nucleations in the 
PNG model of Figure |3J The PNG with sources corresponds to the Ham- 
mersley process with sources and sinks [To] . 

4 Description of the PNG via line ensembles 
4.1 Non- intersecting line ensembles 

The surface height of the PNG model at time T, x i— > h(x,T), does not 
contain anymore the information of the position of the Poisson points, be- 
cause whenever two steps collide, information is lost. Therefore the measure 
induced by the Poisson process on the set of heights is not easy to describe. 
A way of recording the lost information is to extend the model to a multi- 
layer version. This is achieved using the Robinson-Schensted-Knuth (RSK) 
construction [77j . 

Instead of a single line x i— *> h(x, t) evolving in time according to the 
PNG dynamics, one considers a set of lines {x t— > he(x,t)\£ < 0}, with the 
identification ho = h. The initial condition is he(x, 0) = £, x € R, I < 0. 
The deterministic dynamics is identical for all the lines. Only the stochastic 
part differs as follows. The first line follows the PNG dynamics as explained 
in Section [21 If at time t an up- and a down-step annihilate at position x, 
one records the information in the second line in the form of a nucleation 
event at (x,t). This procedure is repeated recursively, that is, the nucleation 
events in the line he-\ correspond to the annihilation events in the line he- 

As for the single line, also for the multilayer version of the PNG model 
there is a geometric point of view. In space-time, when two light-like lines 
generated by the Poisson points meet, they become light-like lines for the 
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£ = £=-1 £=-2 £=-3 



Figure 7: RSK construction up to time t = T. The nucleation events of level 
— 1 are the empty dots and the light-like lines of level 2 are the dotted lines. 
The corresponding line ensemble is represented above. 



second level. More generally, when two light-like lines of level £ meet, they 
continue as lines of level l—l. This is illustrated in Figured By construction, 
the number of lines of level £ crossed along any time-like path from (x, 0) to 
(x, t) is greater than the one of level l — l, for all t < 0. Therefore the set of 
lines has the non-intersecting property 

ht{x,t) >h e -i(x,t) + l (4.1) 

for all x e R, t > 0, I < 0. 
Non-intersecting line ensemble for t = T 

The RSK construction gives us the set of height functions {he, t < 0}. If we 
want to look at the PNG height at fixed time, say t = T, we consider the set 
of height functions x ^ h e (x,T), t < 0. By (J3~H) {h e (-,T),£ < 0} is a set 
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of non- intersecting lines with x h- > Hq{x,T) the surface profile at time T, see 
Figure 

Non-intersecting line ensemble for other space-time cuts 

In some situations, as in the work on the flat PNG [2U| > it can be convenient 
to analyze a line ensemble which corresponds to another space-time cut. 
Consider a continuous and piecewise differentiable path 7 : I — > R x [0, T], 
/ C R an interval. Then the line ensemble corresponding to 7, denoted 
by {Hi,£ < 0}, is given by H e (s) = h^{s)) y s 6 I, £ < 0. It is a non- 
intersecting line ensemble because of (14. lj) . 

4.2 Line ensembles and (real valued) Young tableaux 

We now explain the connection between line ensembles and Young tableaux. 
The height of the PNG surface above x = at time t depends only on the 
Poisson points in the backward light cone of (0, T), i.e., on A( 0) t) = t) G 
Rx R + s.t. \x\ < T — t}. Let us consider the two space-time cuts 71, 72 given 
by 

s^TiW = s), se[0,T], 

s ^ 72 ( s ) = ( S -T,s), se[0,T\. (4.2) 

Denote by < 0}, resp. {Hj, 2 \i < 0}, the line ensemble along 7 1; resp. 

72. From these two line ensembles we construct a pair of Young tableaux 
(Yi, Y2) as follows. Let us start with Yi. Let < si < s 2 < . . . < T be the 
positions of steps in the line ensemble {Hg , I < 0}. Let 4 denote the line 
in which the step at Si occurs. Then has the entry i in row j if the step at 
Si happens in the line Hj^j, i.e., if ^ = 1 — j. Y 2 is obtained in the same way 
with step positions along 72. On the other hand we can define a permutation 
a G Sn, with iV the number of Poisson points in A(o.t), by recording the 
relative positions of the projections of the Poisson points along 71 and 72. 
More precisely, we set (y, z) = (t + x,t — x) and label the Poisson points such 
that Zi < z i+ i. Then, a is the permutation such that y a ^\ < y a u+x\. See 
Figure |H1 for a simple example. At first glance surprisingly, we have 

Yi=7V), Y 2 = Q(a). (4.3) 

But a closer inspection reveals that the multilayer PNG dynamics is just a 
translation of the algorithm in Section l3~3l In particular, this implies that 

H ( l ] (T) = Ef (T) = Ai_* + i (4.4) 
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V{a) 



a = 



(2,3,1) 




Figure 8: Young tableaux and line ensembles. 



where A , Ai, . . . are the length of the rows of the Young tableaux. 

For the PNG droplet studied in [37] , the measure on the Young tableaux 
is the Poissonized Plancherel measure. For flat PNG the measure induced by 
the Poisson points can still be studied if we introduce the symmetric images 
of the Poisson points with respect to the t — axis, see (2Hj for details. 
In [33] a version of the PNG droplet is studied where the Poisson points are 
symmetric with respect to x — 0. These two symmetries correspond to two 
different involutions on the corresponding permutations, denoted by S and 
in 

Length of the rows of Young tableaux. Given the permutation a as 
above, the interpretation of Aj(<r), i = 1, . . . , k, follows from a theorem of 
Greene [HI]. Let, for i < k, a^a) be the length of the longest subsequence of 
cr consisting of i disjoint increasing subsequences. Greene proves that 



In terms of directed polymers, is the maximal sum of the lengths of k non- 
intersecting directed polymers from (0,0) to (t, t), where non-intersecting 
means without common Poisson points. 

Real-valued Young tableaux. From a different point of view, the line en- 
sembles can be regarded as a generalization of Young tableaux for real valued 
entries. Consider a configuration of Poisson points with iV points in A(o,r) 
ordered as above. Then one can run the Robinson-Schensted algorithm with 
the replacements, 



for i — 1, . . . , N. In this way one obtains a pair of Young tableaux with Zi 
as entries. This pair of real-valued Young tableaux contains exactly the same 




(4.5) 




(4.6) 
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information as the line ensembles (H\ , H\ ). More precisely, (resp. Zi) is 
in row j of V (resp. Q) if there is a step in the line H^}- at ?/j (resp. in the 
line h[_a at 2$). 

4.3 Point process associated with the line ensemble 

Our goal is to analyze the surface height h(x,T). Therefore it is natural to 
consider the line ensemble for t = T, {x t— > hi(x,T),£ < 0}. Let us define 
the extended point process j^onRxZ given by 

/ ,\ f 1, if there is a line passing through (x, j), _ x 
= { 0, otherwise. ^ 

r\ T is an extended point process in the sense that, for any fixed x e R, 
t}t(x, •) is a point process on Z j^O]- One can interpret the lines as trajecto- 
ries of fermions in imaginary time (reflecting the non-intersecting condition) 
where the particles, thus the associated point process, evolve following a 
given (stochastic) dynamics [HZj- We now explain the structure of the point 
process r/y and the edge scaling for the PNG droplet and flat PNG. 

Point process for PNG droplet: determinantal 

For the PNG droplet it is enough to consider the line ensemble for x G 
[— T, T], since for x < —T and x > T the lines are straight with hi(±T) = £, 
i < 0. As shown in jSTj, the measure on the configurations is the uniform 
one. Using the transfer matrix method it is shown that t]t is an extended 
determinantal point process, which means that the n-point correlation func- 
tions p( n '(xi,ji] . . . ; x n ,j n ) can be expressed as a n x n determinant. For the 
PNG droplet the kernel is the extended Bessel kernel B T (xi, ji, x 2 , jz), thus 

P {n) (xi, ji] ■ ■ • ;xn,jn) = det (B T {x k ,j k \x h ji)) x < kl < n . (4.8) 

To analyze the surface one defines the edge scaling of the point process, 
compare with (j2.5|) . by 



r^ e (r, s ) = T 1/3 t] T (tT 2/3 , \_2T^l - t 2 T^ + sT 1/3 \). (4.9) 

In the T — > oo limit, this extended point process converges to the extended 
determinantal point process onExR with extended Airy kernel 

f J° oo dAe A ^-^)Ai( Sl -A)Ai( S2 -A), for r 2 > r u 

A(T2, S2\ Ti, Si) = < 

[ - J °° dA e x ^-i) Ai( 5l - A) Ai(s 2 - A), for r 2 < n . 

(4.10) 
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where Ai is the Airy function as given in TJ. In particular for n = r 2 f!4.10|) 
is equal to 

A{ ) = AM^M^MM?A t ( 4.ii) 

s 2 - Si 

the classical Airy kernel [55] . 

Let /i!j? sc be the rescaled height as defined in (j'2.5|) . The joint-distribution 
of /i^? sc at positions r 1; . . . , r n can be written as a Fredholm determinant and 
in the limit of large growth time T 

n 

hm P( fl^fa) < ^}) = det(l - / 1/2 A/ 1 / 2 ) L2((Tl ,..., T „ )xR) (4.12) 
fc=i 

with /j(s) = 0(s — Sj). In this way in |2Z] the convergence of h T ^ sc to the 
Airy process is determined in the sense of finite-dimensional distributions. 

The Airy process 

The Airy process is defined by its finite-dimensional distribution [23 l4*U] . 
For given s\, . . . , s n G R and n < ... < t„ e R, we define / on A = 
{n, ...,r n }xRby f(sj,x) = X( S:j ,oc)(x). Then 

POA(ri) <«!,.. . ,^(r n ) < s n ) = det(l - / 1/2 A/ 1 / 2 ) L 2 (Aid n I . ) 

with A the integral operator with extended Airy kernel. 

The Airy process A was first introduced in [SZ] in the context of the PNG 
droplet. There it was shown that A(t) has a version which is almost surely 
continuous, stationary in t, and invariant under time- reversal. Its single time 
distribution is given by the GUE Tracy- Widom distribution. In particular, 
for fixed t, 

P(A(t) > y) ~ e - y3/H/s for y -> oo, 
P(A(t) < y) ~ e" ls/|3/12 for y -> -oo, (4.13) 
Thus the Airy process is localized. Define the function g by 

Var(^(t)-^(0)) = ^(t). (4.14) 

From jSTj we know that g grows linearly for small t and that the Airy process 
has long range correlations: 

(2t + 0(t 2 ) for \t\ small, 

^ l J \ </(oo) - 2r 2 + <9(r 4 ) for \t\ large. 1 J 

with g(oo) = 1.6264.... The coefficient 2 in front of t~ 2 is determined 
in [2 EH]- The Airy process has been recently investigated and a set of 
PDE's [2] and ODE's [Z21 ESj describing it are determined. 
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Point process for fiat PNG: PfafRan 

The case of the fiat PNG is more difficult and, up to now, the analysis shown 
above is carried out only at fixed position. Because of spatial translation 
invariance, we can choose the reference position to be x = 0. We define the 
edge scaling of the point process as 

v e ^ e {s) = 2- 2 / 3 T 1 /V(0, |_2T + s2 _2 / 3 T 1 / 3 J). (4.16) 

The factor 2~ 2 / 3 is only for convenience, compare with (|2.7j) . It will make 
the connection with GOE random matrices more transparent. 

The structure of this point process is not determinantal as for the PNG 
droplet. In [20J we analyze it and prove that it converges weakly to a Pfaffian 
point process denoted by r/ GOE . More precisely, it is shown that for smooth 
test functions of compact support fi, ■ ■ ■ , f m , m G N, 

m m 

^{U^Vk)) =E(n^ G0E (A))- (4-17) 

k=l k=l 

This point process appears also in the edge scaling limit of the eigenvalues 
of GOE random matrices (see Section 15. 3|) . hence the name ^ GOE . ^ GOE 
is a Pfaffian point process with kernel K GOE , i.e., the n-point correlation 
functions p ( - n \si, . . . , s n ) are given by 

pW( ai> ...,*„)= Pf (K GO *(s h s,))^^ (4.18) 

where K GOE is a 2 x 2 matrix kernel. For an antisymmetric matrix A of even 
dimensions, one has the identity Pf(A) 2 = det(A). The entries of K GOE are 
given by 

POD 

K G ° E ( Sl ,s 2 ) = / d\Ai(s 1 + \)Ai'(s 2 + X)-(s 1 ^s 2 ), (4.19) 
Jo 

/OO -i POO 

dAAi( Sl + A)Ai(s 2 + A) + -Ai( Sl ) / dA Ai(s 2 - A), 

K GOE ( Sl ,s 2 ) = -K GGE (s 2 , Sl ) 

/*oo poo 

^2 G 2° E (si, s 2 ) = - J dA J d//Ai( Sl - A) Ai(s 2 - n) - ( Sl <-> s 2 ), 
with (si <-> s 2 ) standing for the previous term with si, s 2 interchanged. 

5 Random Matrices 

In this section we introduce the Gaussian ensembles of random matrices, the 
Tracy- Widom distributions, and the multi-matrix models. The literature 
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on random matrix is large, the standard reference book is jlTj, see also the 
review [23]. 

5.1 Gaussian ensembles of random matrices 

Let us define the two random matrix ensembles which are linked to the PNG 
droplet and the flat PNG. 

Gaussian unitary ensemble (GUE) 

One defines a Gaussian measure on the set Q of N x N complex Hermitian 
matrices which is invariant under unitary transformations. Let H G Q, then 

P(H)dH = -exp(-Tr(H 2 )/2N)dH (5.1) 

Zj 

where dH = Ylf=i &Hi,i Yli<i<j<N dRe(iJjj)dIm(ifjj) is the product measure 
on the independent coefficients of N. Here and below Z stands for the proper 
normalization constant. The scaling factor 1/2N in (j5.1j) is chosen in order 
to simplify the comparison with the results on the PNG. 

One interesting quantity of random matrices is the distribution of the 
eigenvalues Ai,...,Ajv, which is obtained by integrating over the unitary 
group. The result is 

1 N 

P 2) at(Ai, . . . , X N )dX, • • • d\ N = — |Aat(A)| 2 J] e^^dXj, (5.2) 

3=1 

with A^v(A) the Vandermonde determinant 

A N (X) = det(\l% =1 = J] (A, -A,). (5.3) 

l<i<j<N 

Gaussian orthogonal ensemble (GOE) 

In this case the matrices are real symmetric and the measure is invariant 
under orthogonal transformations. The probability distribution takes the 
same form as in J-LJJ)- with the reference measure dH = rii<t<j<iv 
The integration over the orthogonal group leads to the distribution of the 
eigenvalues Ai, . . . , Ajv 

1 N 

P 1|iV (Ai, • • . , A JV )dA 1 • • • d\ N = — |Ajv(A)| J] e~ x V 2N dX r (5.4) 

3=1 
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5.2 Edge scaling, Tracy- Widom distributions 

One can focus on the largest eigenvalue's distribution when the size of the 
matrices iV is large. The largest eigenvalue, A max , takes value close to 2N and 
its fluctuations are of order N 1 ^ 3 [IH1 EH lEHj • Tracy and Widom determined 
the scaling function of A max with the following result jMUZOI, see also their 
review paper jTT]. Let Fp jN (t) = P/3,Ar(A max < t), then Fp(s) defined by 

F p (a) = Jfin F PtN (2N + sN 1 / 3 ) (5.5) 

exists for (3 = 1, 2. F\ and F 2 are called the GOE and GUE Tracy- Widom 
distribution respectively. More precisely, 

F 2 {s) = exp J (x-s)q 2 (x)dx^ , F 1 (s) = exp ^— ^ J q(x)dx^ F 2 (s) 

(5-6) 

where q is the unique solution of the Painleve II equation q" = sq + 2q 3 
satisfying the asymptotic condition q(s) ~ Ai(s) for s — > oo. 

5.3 Eigenvalues' point process and its edge scaling 

We now define the point process of the eigenvalues Ai, . . . , Ajv- 

GUE point process: determinantal 

We denote by Cat UE the point process on R of the GUE eigenvalues Ai, . . . , Ajv, 



1/2 



i.e., 

N 



C^ UE (A) = ^5(A-A,), AgR, (5.7) 
i=i 

with the abuse of notation 5(X — Xj) for the Dirac measure 5xi- C/v UE is a 
determinantal point process, see for example Chapter 5 of jUj, with kernel 
given by the Hermite kernel 

k h, x y j = PN{x)p N - 1 {y)-p N - 1 {x)p N {y) c _ [x 2 W y AN 
N ' x-y 

where p k (x) = (2nNy 1/4 (2 k k\)~ 1/2 pf(x/V2N) with the (standard) Hermite 
polynomials p^ (x) = e x ' 2 ^e~ x2 . 

The edge scaling of the point process Ctv UE corresponding to ()5.5j) is 

^UE (0 = N l/3 C GUE (2N + eiV l/ 3) _ (5 Q) 
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The limit point process 

„GUE _ ^GUE 



lim V ^ = tj (5.10) 

TV— >oc 



is well defined. It is the determinant al point process with Airy kernel (j4.11jl . 
The GUE Tracy- Widom distribution i*2(s) can also be written as a Fredholm 
determinant 

F 2 (s) = det(l - A) L 2 {(s ^ dx) (5.11) 

with A the integral operator with Airy kernel. For more information on 
determinantal point processes, see for example 



GOE point process: PfafRan 

In the same way as for GUE we define the point process of the GOE eigen- 
values Ctv° E - For even N, Ctv° E is a Pfaffian point process, see for example 
Chapter 6 of For more recent developments on Pfaffian point processes 
see jSHl HH1 • The edge rescaled point process is then given by 

^OE (e) = N l/3 C GOV( 2N + ejV l/ 3) _ (512) 

In the N ^ oo limit, the point process r^ OE converges [7U] to the point pro- 
cess r] GOE which has correlation functions given by ()4.18|) . One consequence 
is that the distribution of the largest eigenvalue can be written as a Fredholm 
Pfaffian (see Chapter 8 of [59J) or Fredholm determinant on the measurable 
space ((s, oo), dx) 

F 1 (s) = Pf ( J - K GOE ) = Vdet(l + JK GOE ) (5. 13) 

with J(x, y) = 5 X)V ^ ^ J j . An interpretation of the Fredholm determi- 
nant as the one on the operator with integral kernel K GOE needs some care 
as pointed out by Tracy and Widom in [Z3j, where they actually prove that 
the finite- N Fredholm determinant converges to Fi(s). 

The GOE kernel being a 2 x 2 matrix is not uniquely defined. In fact 
using Pi(X t KX) = det(X) Pf (K) [HE] one obtains a family of kernels K 
leading to the same correlation functions. For example, the kernels reported 
in [2H1ISS] differ slightly from ()4.19|1 . but they are equivalent since they yield 
the same point process. 

Very recently a new formula for the GOE Tracy- Widom distribution Fi 
was established. Let B(s) be the operator with kernel B(s)(x,y) = Ai(x + 
y + s). It then holds 

Fi(s) = det(l -B(s)) (5.14) 

where the Fredholm determinant is on L 2 (R + ). This is proven in [221 starting 
from the work on the flat KPZ growth [63J. 
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5.4 Multimatrix model: Dyson's Brownian motion 

Dyson noticed that the distribution of eigenvalues (|5.2J) and ()5.4j) is 
identical to the equilibrium probability distribution of the positions of N 
point charges, free to move in R under the forces deriving from the potential 
U at inverse temperature f3, with 

1 N 

U{xx,...,X N )=- ln \ X i~ X i\ + l^n^2 X l ( 5 - 15 ) 

In the attempt to interpret the Coulomb g dynamical system Dyson 

considered the positions of the particles in Brownian motion subjected to 
the interaction forces — V£7 and a frictional force / (which fixes the rate 
of diffusion, or equivalently, the time scale). He showed that in terms of 
random matrices, this is equivalent to the evolution of the eigenvalues when 
the independent coefficients of the matrix H = H{t) evolve as independent 
Ornstein-Uhlenbeck processes given by 

P(H(t) = H\H(0) = H )dH = i exp ) (5.16) 

with q = exp(—t/2N). The evolution of the eigenvalues satisfies the set of 
stochastic differential equation 

<%« = ( - Imo + f t jm^m)* + *>(«) , i = i, .... iv, 

(5.17) 

with {&j(t), J = 1, ■■■,N} a collection of iV independent standard Brownian 
motions. The parameter f3 is 1 for GOE and 2 for GUE. We refer to the sta- 
tionary process of \5.1 r /\j as Dyson's Brownian motion (for the eigenvalues). 
Note that for (3 > 1 the process is well defined because there is no crossing 
of the eigenvalues, as proved by Rogers and Shi [BTj . 

5.5 GUE extended point process 

The GUE eigenvalues {A,(t), % — 1, . . . , iV} define an extended point process 
on R + x R, 

v 

C^MH^A-A^)). (5.18) 
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It is a determinantal point process JH] with kernel given by the extended 
Hermite kernel 1 , 



K^(t 2 , X 2 ; ti, Xi) 



E fc =- W e fe ^^)/ 2 V(AiK(A 2 )e-( A ? +A i)/^, t 2 > t u 

-Er=oe fc( ' 2 -' l)/2 V(AiK(A 2 )e-^^)/^, t 2 < t lt 

(5.19) 



where pk(x) = p^ +k (x/ '^/2N)(^/2TrN2 k k\)~ 1 ^ 2 and pf(x) the standard Her- 
mite polynomials. 

The GUE Dyson's Brownian motion is stationary and the edge scaling is 
the following. Let U = 2r;iV 2 / 3 , \ = 2N + SiN 1 / 3 , i = 1,2. In the limit of 
large N, converges in the edge scaling limit to the extended Airy kernel 

lim N 1/3 K^(2t 2 N 2/3 ,2N + s 2 N 1/3 -2t 1 N 2/3 ,2N + s 1 N 1/3 ) = A(r 2 , s 2 ; n, 

N^oo 

(5.20) 

This convergence is uniform for s%, s 2 G [a, 00) for any fixed a £ R, see for 
example Appendix A. 7 of 21. Moreover, since it has (super-)exponential 
decay for large s\,s 2 — > 00, one deduces that the largest eigenvalue converges 
to the Airy process in the sense of finite-dimensional distributions, 

lim A^- 1 / 3 (A max (2riV 2 / 3 ) - 2N) = A{r). (5.21) 



6 PNG model, Young tableaux, and random 
matrices 

6.1 PNG model 

Although the PNG model and the random matrix ensembles describe systems 
completely different, we can associate some point process via the multilayer 
extension and the integration over the symmetry groups, respectively. The 
point processes happen to have the same mathematical structure: 

(a) determinantal point processes for PNG droplet and GUE random ma- 
trices, 

(b) Pfaffian point processes for flat PNG and GOE random matrices. 

Moreover we have seen that, in the appropriate scaling limit, they converge 
to the same limit objects. Thus the distribution of the largest eigenvalue and 

1 According to |73] this was already described in a MSRI lecture (2002) by Kurt Jo- 
hansson. It can be derived for example using Theorem 1.7 of |39| . details can be found for 
example in Appendix A. 6 of |21| . 
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the height of the PNG at fixed position have the same fluctuations in the 
asymptotic limit. 

One natural question is whether such a connection still exists at the level 
of joint-distribution (extended point processes). The extended point process 
for the PNG model is naturally defined from the line ensemble. As explained 
in Section 15.41 one defines the so-called Dyson's Brownian motion (multi- 
matrix ensembles) and associates to it an extended point process. For GUE 
random matrices it turns out that the extended point process has, in the 
limit of large matrix dimension N, the extended Airy kernel, the same as for 
the PNG droplet. As a consequence, the (droplet) surface height and the 
evolution of the largest (GUE) eigenvalue are described by the same process 
in the asymptotic limit [22] : the Airy process. 

For GOE random matrices the answer is not known. It is conjectured that 
the surface height of the flat PNG and the evolution of the largest (GOE) 
eigenvalue converge to the same limit process. The result of Ferrari [20] makes 
this conjecture more plausible. In fact, we now know that, not only h(0,T) 
in the limit T — > oo and properly rescaled is GOE Tracy- Widom distributed, 
but also that the complete point process r\ T converges to the edge scaling 
of Dyson's Brownian motion with (3 = 1 for fixed time. For (3 = 1 Dyson's 
Brownian motion this structure has not been unraveled and the question 
remains open. Very recently, Sasamoto [HB] obtained a process in the context 
of the totally asymmetric exclusion process, which, by universality, should 
describe also the surface height of the flat PNG. 

From the point of view of multi-matrices, the difficulty is the fact that 
the Harish-Chandra/Itzykson-Zuber formula doesn't have a particularly nice 
analogue for symmetric matrices. From the point of view of the SDE (I5.17J1 
the problem lies in the different factor in front of the drift term, which seems 
to make it impossible to regard the process as Doob transforms of N inde- 
pendent processes pjjj . 

There are other models with the same mathematical structure and show- 
ing the same limit behaviour (and scaling functions): the 3D-Ising cor- 
ner [5*2*1 |2I] and the Aztec diamond jHHl E3 which are linked via domino 
tilings, the vicious random walks and the non-colliding Brownian particles, 
see for example 0H1E2]; and the totally asymmetric exclusion process [2~""*. l63j. 
Some of the problems presented above can be solved by using orthogonal 
polynomials point of view instead of the more geometric line ensembles, see 
for example the review [13]. In this context the relevant processes are the 
Schur process [!52] and its Pfaffian analogue [3*4"| IT*"]. 
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6.2 Young tableaux 

In Section 14.41 we described the connection between the length of the row of 
the Young tableaux and the point process of the multilayer PNG. Above we 
explained the connection between the multilayer PNG and the eigenvalues 
of random matrices. Putting the two arguments together it follows that 
there is a connection between Young tableaux and random matrices. After 
having seen that the length of the first row of Young tableaux under the 
Plancherel measure (jH.5j) and the largest eigenvalue have the same statistics 
in the asymptotic limit, it was conjectured by Baik, Deift, and Johansson that 
the connection would extend to the top rows, respectively top eigenvalues. 
They proved it for the second row in [5] and the first proof for all top rows 
is due to Okounkov [ST]. Other proofs are in |3T[ ITS] . In the work on 
flat PNG [20] the point process with symmetric images was studied, which 
corresponds to the involution S for permutations. In this case, the measure 
on the Young tableaux becomes Pl^(A) = d\, A G Yjp , being the 
restriction of Yjv to the Young tableaux with Aj even for all i. From the 
proof for the flat PNG [20] it follows that the top rows of Yjp and the top 
eigenvalues of GOE have the same limit joint-statistics. 

We first learned of the connection between the Young tableaux and the 
random matrix ensembles for the GOE and GSE random matrix ensembles 
by Sasamoto [02]; where he deduced the connection starting from the PNG 
half-droplet [33J. For the GOE case the restriction that all A« have to be even 
is not required. This result is also obtained by Forrester, Nagao, and Rains 
in [20] using an approach with orthogonal polynomials. 
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